Orthogonality catastrophe in a one— dimensional system of correlated electrons 
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We present a detailed numerical study of the orthogonality catastrophe exponent for a one- 
dimensional lattice model of spinless fermions with nearest neighbor interaction using the density 
matrix remormalization group algorithm. Keeping up to 1200 states per block we achieve a very 
great accuracy for the overlap which is needed to extract the orthogonality exponent reliably. We 
discuss the behavior of the exponent for three different kinds of a localized impurity. For comparison 
we also discuss the non-interacting case. In the weak impurity limit our results for the overlap con- 
firm scaling behavior expected from perturbation theory and renormalization group calculations. 
In particular we find that a weak backward scattering component of the orthogonality exponent 
scales to zero for attractive interaction. In the strong impurity limit and for repulsive interaction 
we demonstrate that the orthogonality exponent cannot be extracted from the overlap for systems 
with up to 100 sites, due to finite size effects. This is in contradiction to an earlier interpretation 
given by Qin et al. based on numerical data for much smaller system sizes. Neverthless we find 
indirect evidence that the backward scattering contribution to the exponent scales to 1/16 based on 
predictions of boundary conformal field theory. 

PACS numbers; 71.10.Pm, 72.15Nj, 79.60.-i, 72.15.Qm 



I. INTRODUCTION 

Anderson's "orthogonality catastrophe" is an impor- 
tant and universal piece of metallic physics with far 
reaching consequences for the response functions of 
many-electron systems. It is well known that for non- 
interacting electrons the overlap between the many-body 
ground states with and without a local impurity poten- 
tial vanishes in the thermodynamic limit [1]. For a large 
but finite system of length L the overlap tends to zero 
as a power law L^", defining an orthogonality exponent 
(OE) a. This orthogonality catastrophe is closely related 
to the x-ray edge problem in metals [2,3] , and in this con- 
text the exponent a can often be measured directly. 

As was first observed by Anderson [1] the OE for a 
spherically symmetric impurity in a free electron gas 
(in arbitrary dimension) is uniquely determined by the 
square of the scattering phase shifts of the infinite sys- 
tem at the Fermi energy Ep. The proof of this result 
for more general potentials is far from trivial [2,4,5]. For 
spinless one-dimensional (ID) Fermions we can write 



with the phase lpt of the transmission coefficient and the 
absolute value \R\ of the reflection coefficient taken at 
the Fermi energy. If we introduce new phase shifts (5y /{, 
given by 
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we can identify the forward scattering phase shift <5/ 
with and the backward scattering phase shift E^, with 
arcsin Expressed in terms of (5j/f, or and (/?t, 
the OE is given by 
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with the phase shifts of the even and odd scattering chan- 
nels. The remaining problem is the calculation of the 
phase shifts for a given impurity potential. In the follow- 
ing 8 always denotes the phase shift at the Fermi energy. 

In one dimension the phase shifts can generally be ex- 
pressed in terms of the reflection and transmission coef- 
flcients of an incident wave packet [6] 



deflning the forward (a/) and backward (ah) scattering 
contributions to the OE. 

In the lowest order perturbation theory in the impurity 
strength (the Born approximation) it is always possible 
to express the scattering phase shifts 5^/^ i^i terms of the 
matrix elements of the impurity term W in the Hamil- 
tonian. If we denote the one-particle eigenstates of the 
impurity free Hamiltonian with energy Ep by \^kp), we 
obtain 
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Sejo = t,VptT arcsin{|i?|}] , 
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Correlation effects are very important in one dimen- 
sion, where it is known that even an infinitesimally small 
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interaction between electrons changes their low energy 
properties from those of a Fermi to those of a Luttinger 
liquid (LL) [7]. The low energy excitations of a LL are 
collective and bosonic, so the concept of phase shifts in 
fermionic one-(quasi-)particle wave functions is clearly 
no longer applicable, contrary to the situation in interact- 
ing Fermi liquids where this is still possible. It is therefore 
interesting to ask whether the orthogonality catastrophe 
in a ID system is also dramatically altered by electron- 
electron interaction. 

In this article we analyze the problem by numerically 
calculating the overlap 

0=\{Eo\E^)\ (7) 

and the difference in the ground state energies 

AE=\E^-Eo\ (8) 

between the ground states of the system with {\Eq)) 
and without (|ii^o)) the impurity for three different kinds 
of impurity in a lattice model of spinless interacting 
Fermions. We distinguish between the weak and strong 
impurity limit. The results are compared with the pre- 
dictions of perturbation theory and the application of 
renormalization group (RG) ideas. Large system sizes 
and very accurate calculations are found to be necessary 
to reliably determine a. Th(^s(^ can be achieved by us- 
ing the density matrix renormalization group algorithm 
(DMRG) [8,9]. This approach makes it possible to calcu- 
late the ground state properties of an interacting system 
with an accuracy comparable to that of exact diagonal- 
ization, but for much larger system sizes. Using this tech- 
nique, we are able to work with chains of up to 100 sites, 
keeping up to 1200 states per block and performing 7 
finite lattice sweeps in the DMRG procedure. 

The behavior of ID correlated electrons in the presence 
of a localized impurity is of general interest because the 
physics of the problem manifests itself in a number of dif- 
ferent contexts, e. g. the conductivity of quasi ID "metal- 
lic" materials [10], tunneling through a constriction in the 
fractional quantum Hall regime [11] and Friedel oscilla- 
tions in the charge density of highly anisotropic systems 
[12 14]. Given the rapid progress in the fabrication of 
quasi ID metallic systems (highly anisotropic materials 
or quantum wires) it should in the future be possible to 
measure the x ray response of these systems directly, and 
so fit a to experiments. 

Considerable progress has been made in the under- 
standing of the intcTplay between impurity and interac- 
tion in one dimension since the seminal work of Kane 
and Fisher [10]. Kane and Fisher discussed the prob- 
lem of a single impurity in a ID wire within the effective 
low energy field theory of ID correlated electrons, using a 
perturbative RG approach. If one integrates out high mo- 
mentum degrees of freedom any weak 2k p component of 
the impurity potential increases for repulsive interaction 
and scales to zero for attractive ones. Thus backward 



scattering of electrons is relevant for repulsive interac- 
tion and irrelevant for attractive interaction. Kane and 
Fisher also discussed the dual problem of a weak hopping 
between the two ends of an open chain. The strength of 
this weak link flows to zero for repulsive interaction, but 
is relevant for attractive interaction. 

It is tempting to conclude that, for repulsive interac- 
tion, the (weak) 2kF component of the impurity potential 
scales to infinity and therefore the open chain fixed point 
(open chain interpretation). However it is not immedi- 
ately obvious that a single impurity in a LL really does 
correspond to a "cut wire" . 

In order to shed more light on this issue Kane and 
Fisher calculated the conductivity of a ID wire with a 
localized impurity exactly for a special value of the LL in- 
teraction parameter. Their result is consistent with this 
interpretation of the RG results [10]. Later authors calcu- 
lated the conductivity for arbitrary interaction strength 
[15,16], again obtaining results consistent with the open 
chain interpretation. Further supporting evidence for 
this picture was given by a numerical calculation of the 
conductivity [11]. Similar arguments hold in the dual 
situation of (strong) hopping between two semi-infinite 
LL's, attractive interaction, and the periodic chain fixed 
point (periodic chain interpretation). 

However the question of how the RG results should be 
interpreted in the context of the OE remains controversial 
[17-22]. 

A number of attempts have been made to calculate 
the orthogonality exponent (or equivalently the x-ray 
edge exponent) of a ID interacting electron gas with 
strong impurities in the light of Kane and Fisher's RG 
results [17-19]. All of these approaches are based on the 
Tomonaga Luttinger (TL) Hamiltonian, the effective low 
energy continuum field theory for most models of ID cor- 
related electrons in the absence of impurities. The TL 
Hamiltonian is written in terms of distinct left and right 
moving fermion fields 'ip±{x). The low energy physics 
of a given microscopic model can be mapped onto a TL 
model by determining the LL parameter K and the renor- 
malized velocity v of the model in question, in terms of 
its microscopic parameters [23]. Impurity scattering is 
given by two types of term in the Hamiltonian. Forward 
scattering can be written linearly in the density of right 
and left moving Fermions. If only this kind of scatter- 
ing is present then the calculation of the overlap within 
the field theoretical model is straightforward (see Section 
lie). Backward scattering couples right and left moving 
fields and results in a complicated sine-Gordon like term 
in the field theory [10]. 

To calculate the overlap within the field theory in the 
presence of backward scattering further approximations 
are needed. To solve the problem for repulsive electron- 
electron interaction Gogolin [17] and Prokof'ev [18] used 
a strong backward scattering approximation, replacing 
the cosine of the boson field by a term quadratic in the 
field. This replacement is consistent with the open chain 
interpretation of the RG results and leads to a backward 
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scattering contribution 



to the OE, which is completely independent of the 
strength of the interaction or the size of the bare back- 
ward scattering potential. 

This result is frequently interpreted in terms of phase 
shifts. According to the open chain interpretation of 
Kane and Fisher's scaling analysis outlined above, the 
2kF component of the impurity is renormalized to in- 
finity. This leads to perfect reflection (reflection coeffi- 
cient R = 1) and a backward scattering phase shift of 
7r/2. Introducing this phase shift into Eq. (4) gives Eq. 
(9). As there arc no fcrmionic single particle scattering 
wave functions which are eigenfunctions of the interact- 
ing Hamiltonian, this interpretation only has a mean- 
ing if one extends the concept of phase shifts. We can 
do this by regarding the phase shift at the Fermi en- 
ergy as the number determining the boundary condition 
between in- and outgoing field operators according to 
V'out = exp (2z(5)?/'in- In this sense the above interpre- 
tation does still have a physical meaning [24] . 

Beside the direct calculation of the overlap, an alterna- 
tive method of calculating the OE is given by boundary 
conformal field theory (BCFT) [24,25]. BCFT relates the 
properties of the boundary changing operator associated 
with the introduction of a single impurity in a (semi-) 
infinite ID interacting Fermi system to the finite size 
spectrum of an equivalent field theory on a line of length 
L. For non-interacting electrons and for interacting elec- 
trons with only a forward scattering impurity component 
the difference AE of the ground state energies of the fi- 
nite size system with and without the impurity can be 
calculated directly from a mode expansion in the effec- 
tive low energy field theory on the line. In the absence 
of bound state effects AE is given by 

lim LAE= lim LAe + 2'Kva, (10) 

L—^oo L—^oo 

where Ae is a bulk constant and v the renormalized 
charge velocity of the field. For interacting electrons this 
result enables us to calculate a/ directly from numeri- 
cally determined ground state energies [24,26]. 

An alternative derivation of Eq. (1) for non-interacting 
electrons based on Eq. (10) has been given [25]. This ap- 
proach has also been used in conjunction with the Bethe 
Ansatz to obtain predictions for the orthogonality expo- 
nent in lattice systems with integrable defects [27]. Ad- 
ditionally it is possible to calculate the ground state en- 
ergy difference AE in the presence of an infinitely strong 
backward scattering potential. If one assumes that the 
scaling of the 2fci? component of the potential for re- 
pulsive electron-electron interaction finally leads to an 
open boundary condition fixed point one again obtains 
ttb = 1/16. 

It is important to notice that all the theoretical results 
discussed so far are based on, or are equivalent to, the 



open chain interpretation of the results of the perturba- 

tive RG. 

Recently the overlap O Eq. (7) has been discussed 
within the low energy continuum field theory for the spe- 
cial LL parameter K = 1/2 corresponding to a strong 
repulsive interaction [28,29]. At this point it is possible 
to calculate the overlap by bosonization and refermion- 
ization. The result is again consistent with at = 1/16. 

In Refs. [20] and [21] Orcg and Finkel'stein discuss the 
tunneling density of states near a localized impurity and 
the x-ray edge exponent based on the mapping of the TL 
model onto a Coulomb gas model. They obtain results in 
contradiction with results obtained by all other authors. 
Whether their calculations suffer from a flaw in the an- 
ticommutation relations of the Fermion field operators is 
still a matter of contention [22,28-31]. 

In Refs. [30] and [31] Qin et al. use the DMRG method 
to calculate the overlap for a very special type of localized 
impurity. We will later comment on the relation between 
their results and ours. 

The remainder of this article is organized as follows: 
In Sec. II we introduce the model considered, discuss an- 
alytical results for the OE in the non-interacting case, 
and derive the results for the OE in lowest order per- 
turbation theory in the impurity strength, based on the 
effective low energy field theory for interacting electrons. 
In order to obtain the OE it is necessary to extrapolate 
finite size data for the overlap and the energy difference 
to the thermodynamic limit. This leads to an extrapola- 
tion error. An additional error is introduced because the 
DMRG gives only an approximate result for the ground 
state properties. We discuss both errors carefully in Sec. 
Ill by calculating the overlap O and AE numerically for 
the non interacting model and comparing with the exact 
results derived in Sec. II B. 

For non-interacting electrons a straightforward numer- 
ical calculation of O and AE is possible for very large 
system sizes, if we take the Slater determinant structure 
of the ground states into account. The comparison of ex- 
act results for non-interacting electrons makes it possible 
to access the errors involved in the DMRG procedure. It 
turns out that this kind of error analysis is crucial for 
the interpretation of our results. In Sec. IV we present 
our numerical results and analyze the data with respect 
to the expected scaling behavior. Finally we summarize 
our results in Sec. V. In the Appendix we briefly describe 
how to obtain perturbative results in the weak impurity 
and weak hopping limits. 

II. MODELS AND ANALYTICAL RESULTS 

A. The Hamiltonian 

In our investigation we consider a model of N spinless 
Fermions on a ID lattice of M sites with site dependent 
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on-site energies Si , nearest neighbor hopping t and inter- 
action V, described by the Hamiltonian 



i=l i=l 
M 

+V ^ n.m+i + {b-l) -t + CiCm) 



i=l 

+{b' - 1) [VuMni] , 



(11) 



where c) ' denotes the Fermion lowering (raising) opera- 



(t) 

i 

tor at the site i, and the related density operator. In 
all sums we identify M + 1 = 1 . By choosing b = b' = 1 
and Si = for all i we obtain an impurity free Hamil- 
tonian for interacting Fermions with periodic boundary 
conditions (PBC) . It is well known that this Hamiltonian 
displays LL behavior, and for the case of a half filled band 
the LL parameter K and the renormalized velocity v can 
be obtained from the finite size corrections [32] to the 
Bethe Ansatz [33] solution 



V = t 



TTsin (2?7) 



4ry' 

where r] parameterizes the interaction 
V = -2icos(2r/). 



(12) 
(13) 

(14) 



For V > (repulsive interaction), K < 1, and for F < 
(attractive interaction) K > 1. In this article we will fo- 
cus on the half filled case. According to the LL theory of 
ID correlated electrons [7] K and v completely determine 
the low-energy behavior of the Hamiltonian. It is known 
that in the absence of impurities the model described 
by Eq. (11) undergoes a transition into a charge density 
wave ground state at V/\t\ = 2, and into a phase sepa- 
rated state at V/\t\ = —2. We are therefore hmited to 
parameters V/\t\ within (—2, 2) leading to if G (1/2, oo). 

By varying the parameters Si , b and b' we can introduce 
different kinds of localized impurity to the Hamiltonian 
Eq. (11). In this article we consider three kinds of im- 
purity. Setting £i = £, £j = for z 7^ 1 and b = b' = 1 
creates a site impurity on one of the lattice sites. With- 
out loss of generality we have chosen the site 1 as the 
impurity site. Physically one might think of this kind of 
term as an impurity ion sitting on one of the sites of the 
regular lattice, or as the core hole generated by remov- 
ing one (or more) of the core electrons of the atom at 
site 1 (x— ray edge problem). In this article we will focus 
on repulsive impurities (e > 0) to avoid complications 
due to bound states. A hopping impurity is obtained by 
choosing = for all z, b' = 1 and < & < 1 and 
might be created by an interstitial between lattice sites 
M and 1 which changes the hopping amplitude but not 
the interaction between them. 



Qin et al. [30,31] recently discussed a third type of im- 
purity. These authors worked with the Heisenbcrg XXZ 
spin 1/2 chain, which can be mapped onto a model of 
spinless Fermions by a Wigner-Jordan transformation 
[34]. In terms of the Heisenberg model, their impurity 
is generated by weakening the coupling of the x-, y- and 
z-components of the spin between sites M and 1 by an 
equal amount. In the Fermion basis, this impurity corre- 
sponds to an exotic combination of modified hopping and 
interaction between the sites M and 1, accompanied by 
site impurities at the sites M and 1. Their model is equiv- 
alent to Eq. (11) with ei = Em = -{b - l)V/2, Si = 
for all the other i, b — b' and < 6 < 1. For comparison 
we will also discuss this type of impurity, and denote it a 
bond impurity. It is equivalent to the hopping impurity 
in the non-interacting limit V = 0. 

The ground state of the Hamiltonian Eq. (11) with- 
out impurities is degenerate for even particle number N. 
There are several ways of overcoming the problems asso- 
ciated with this degeneracy and (physical) even odd ef- 
fects. One straightforward solution is to calculate O and 
AE only for odd N. However, in order to have access to 
more data points we have chosen a different approach in 
the majority of our numerical calculations. For even A'^ 
we altered the hopping between the sites M and 1 from t 
to —t. This lifts the degeneracy of the ground state, and 
prevents even-odd effects irrelevant to the orthogonality 
catastrophe considered in this article. One can imple- 
ment this in the Hamiltonian by the direct substitution 



(6-1) 



- exp 



M 



(15) 



in the second line of Eq. (11). If one works with the 
Heisenberg Hamiltonian and imposes PBC's for the spins 
this additional phase factor emerges as a natural conse- 
quence of the Wigner-Jordan transformation [34] . 



B. The orthogonality exponent for non— interacting 
electrons 

As shown by Eq. (4), we can determine the OE a 
for non-interacting electrons by solving the one-particle 
scattering problem and calculating \R\ and ipT- For 
the types of impurities considered in this article this is 
straightforward and can be done following e. g. Ref. [35]. 
For the site impurity we obtain 



1 

2^ 



arcsm 



V^4i2 + £2 



(16) 



The forward and backward scattering contributions ctf/t 
are both given by one-half of a, as <^t = — arcsin \R\, and 
\R\ is given by the argument of the arcsinc in Eq. (16). 
In preparation for the discussion of the OE in terms of 
an effective low energy field theory in the next section. 
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wc also calculate 5f /jj in the Bom approximation and at 
half filling. This will be necessary in order to relate the 
impurity parameters of the microscopic lattice model to 
the impurity potentials in the field theory. We denote the 
value of the OE obtained from the Born approximation 
as . Following Eqs. (5) and (6) we obtain {SjY = 
(Jf )2 = {e/vpY with the non- interacting Fermi velocity 
Vf = 2|^|. Inserting this result into Eq. (4) leads to the 
same result as an expansion of Eq. (16) to lowest order 
in £ 



a 



- - 

27r2 \vF 



(17) 



The next order corrections to are of the order [e/vpY. 

To simplify the calculation in the case of the hop- 
ping impurity we concentrate on the half filled band with 
Ep = 0. In this case 



1 . 2 f /3 - /3V2 

a = — - arcsm < 

47r2 \i_/3 + /32/2 



(18) 



with (} = (1 — 6). For the hopping impurity at half fill- 
ing the forward scattering contribution vanishes, e. g. we 
obtain 



5f = ipT = 



(19) 



and \R\ is again given by the argument of the arcsine in 
Eq. (18). For a weak impurity (small /3) an expansion 
leads to 



(20) 



and ((5^) = 0^- In contrast to the site impurity the 
next order corrections to the Born approximation are of 
the order 0^ and therefore more important than for the 
site impurity. This observation will be important in what 
follows. As the impurity strength in the lattice model is 
not determined by a local impurity potential, the calcu- 
lation of the backward scattering phase shift in the Born 
approximation following Eq. (6) is a little bit more in- 
volved, but also leads to Eq. (20). 



C. Bosonization and perturbation theory in the 
impurity strength 

In the low energy, long wavelength limit, the impurity 
free Hamiltonian can be written in terms of a bosonized 
continuum field theory [23]. The new effective Hamilto- 
nian is given by 



Ho 



= — < 

277 J_i/2 



dx { K'k''W{x) 



dx 



(21) 



and the canonical conjugate bosonic fields $(.x) and n(.T) 
can be expressed in terms of the densities p± of right and 
left moving Fermions [23] 



n(a;) = p+{x) - p-{x), 
$(x) = -i^5]ie-«-[p+(g)+p_(«)] 



TTx r 



(22) 



(23) 



where the number operators N± measure the particle 
number with respect to a fixed ground state. In the non- 
interacting model K = 1 and the velocity v is given by 
the non-interacting Fermi velocity vp- 

A forward scattering impurity can be expressed as 



L/2 



L/2 



Wf{x)p{x), 



with the total density p{x) = [/5+(.'z;) -I 
local potential Wf{x). Eq. (23) implies 



p{x) 



1 d^x) 

TT dx 



(24) 

P- (x)] and the 
(25) 



Since W/ is linear in the p±{q), the ground states of Tig 
and Ho+Wf are related by a unitary transformation and 
the overlap can be calculated explicitly 



Of = exp < — 



1 K 



E 

9>0 



(26) 



Taking the limit L — > oo we can determine the OE [12,36] 



af 



47^2 



(27) 



Anderson's orthogonality catastrophe is therefore still 
present in an interacting system with forward scattering, 
but the OE is modified. Eq. (27) trivially reproduces the 
non-interacting limit, and we conclude from comparison 
with Eq. (5) that the field theoretical description always 
gives the Born approximation in Wf for the OE [37] . We 
therefore anticipate that for V ^ 0, Eq. (27) also holds 
only for weak forward scattering. 

This raises the general question of how the specific 
impurities (e. g. bond, site, hopping) of a given micro- 
scopic model should be mapped onto the impurity po- 
tential Wf{q) of the effective field theory. Comparing 
Eq. (27) for y = with Eq. (4) we find 



mm 



(28) 



where 6f is the exact non-interacting forward scattering 
phase shift of the underlying microscopic model. Within 
this mapping we obtain from the field theoretical Born 
approximation the exact result for a f of the microscopic 
model. But it is not clear whether this identification 
also gives a systematic improvement in the interacting 
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case. From a systematic point of view wc can only iden- 
tify the impurity parameters in lowest order in perturba- 
tion theory [Wf{Q)/vF]'^ = (Sf)'^, where is the non- 
interacting microscopic forward scattering phase shift in 
the Born approximation. 

For small |W/(0)| we can therefore write 



[10], and is related to the anomalous time dependence 
of the core hole Green's function discussed in Ref. [17]. 
It is not obvious from Eq. (32) that we will recover the 
non-interacting form 



w2 

Ob = 1-ai^ -af ln(L) + ... 



(33) 



af{V)=K^af{V = 0). 



(29) 



Having calculated o:J{V = 0) for our microscopic site 
impurity in Eq. (17), a^{V) is determined by Eq. (29). 
We will later compare the values of a/ predicted by this 
equation with our numerical results for weak site impu- 
rities. 

To the best of our knowledge there are no analytical re- 
sults for the case of strong forward scattering and V ^ 0. 

It is important to note that the inclusion of the forward 
scattering term Eq. (24) in the Hamiltonian breaks its 
particle hole symmetry. As the lattice Hamiltonians Eq. 
(11) for the hopping and bond impurities at half filling are 
particle-hole symmetric, in this special case the mapping 
of the microscopic model onto an effective field theory 
cannot give rise to such a term. In this sense the forward 
scattering contribution to the OE vanishes for interacting 
electrons as well, which simplifies the interpretation of 
some of our results. 

For purposes of comparison with numerical data for 
finite size systems, it is interesting to evaluate Eq. (26) 
with finite L. Provided that we can expand 11^/(9)1^ in a 
power series in q, the Euler summation formula [38] gives 



in the limit K 1. By taking the explicit V dependence 
of the coefficients ai{V) into account, we show in the 
Appendix that this is in fact the case. 

For repulsive interaction {K < 1), the expansion Eq. 
(32) breaks down at large system sizes, even if the bare 
backward scattering potential is weak. The same is true 
in the non-interacting case, but the logarithmic diver- 
gence there is replaced by a faster power law divergence. 
As we do not find the usual behavior Eq. (33) in the inter- 
acting case, we expect that the orthogonality catastrophe 
in the presence of backward scattering and interaction is 
drastically altered. For attractive interaction {K > 1) 
and in the lowest order perturbation theory Ob tends to 
a constant. Provided that this holds also in higher or- 
ders, there is no orthogonality catastrophe, i. e. the two 
ground states are not orthogonal to each other. Where 
the bare |Wfc| is small, wc expect to find the scaling be- 
havior Eq. (32) in our numerical data for the overlap, for 
intermediate system sizes. Terms of higher order in l/L 
will be important for very small systems. 

As in the forward scattering case a systematic mapping 
between the microscopic backward scattering parameter 
and Wb is only possible within perturbation theory. In 
lowest order we can identify [M^6(0)/i;f]'^ = {^b )"^! where 



O = exp In (L)+ao + aiL-i+a2i-' + ...}, (30) is the non -interacting microscopic backward scatter- 



with impurity strength and interaction dependent con- 
stants a,;. In the Appendix we will argue that in the 
presence of backward scattering the finite size corrections 
for F = are also given by integer powers of 1/L. 

The backward scattering contribution to the low en- 
ergy field theory is given by 



L/2 

m= I Wb{x) vi(a;)V-(a;) + h.c. 

-L/2 L 



(31) 



where the field operators (.x) can be expressed in terms 
of the Boson fields n(x) and $(x) [7,23]. In the Appendix 
we derive an expression for the overlap perturbatively 
in the backward scattering potential. If only backward 
scattering is present [12] 



ing phase shift in the Born approximation. 

In Ref. [10] Kane and Fisher discuss the strong im- 
purity limit by introducing a weak hopping in an open 
chain and mapping this problem to a field theory dual 
to Eq. (21). This mapping makes it possible to treat the 
problem perturbatively in the weak hopping. In our case 
this corresponds to a calculation of the overlap between 
the state with the weak hopping (the large hopping or 
bond impurity) and the chain with open boundary con- 
ditions (OBC). To distinguish this overlap from the one 
discussed previously, we denote it 0° , and the related OE 
a°. As discussed in the Appendix we can calculate O^ 
perturbatively in the weak hopping limit in a similar way 
as Ob and the result is given by an expression equivalent 
to Eq. (32) but with K replaced by l/K 



Ob 



l-ai(y)^-a2(y)-^L^-2^ + ..., 



a,{V)'f,-~a2{V)%L^-''^ + ..., (34) 



(32) 



where the ai{V) are interaction dependent constants, Wb 
is a measure of the backward scattering strength and the 
dots denote terms of second order in Wb/vp falling off as 
1/L or faster and terms of higher order in Wb/vp- This 
result is consistent with the RG equations given in Ref. 



where tb is now a measure of the weak hopping. The roles 
of repulsive and attractive interaction are interchanged 
and Ol tends to a constant for repulsive interaction. We 
will verify Eq. (34) numerically in Sec. IV B. 

BCFT predicts that ab between the OBC and the PBC 
ground states is 1/16, and that this result is independent 
of V [24] . Bearing in mind all the results discussed so far, 



6 



it is therefore tempting to conclude that Q;;,(y < 0) = 
and aii(V > 0) = 1/16 independent of the bare impu- 
rity strength. We want to emphasize that these general- 
izations of the above perturbative results are consistent 
with the periodic and open chain interpretation of the 
perturbative RG results, but need to be verified. For our 
models we will discuss this question in Sec. IV. 

In the next section we will discuss the two major 
sources for errors in our numerical calculation of the OE: 
The extrapolation error, and the error due to the approx- 
imative nature of the DMRG. In this section we will also 
introduce the method used to analyze the numerical data 
for O and AE in order to obtain a. 



III. ACCURACY OF THE NUMERICAL 
RESULTS 

A. Extrapolation to infinite system size 

For non-interacting electrons, the ground states with 
and without the impurity are Slater determinants of sin- 
gle electron states, and the overlap can be written as the 
determinant of a matrix whose elements are given by the 
overlaps between these single electron wave functions [1]. 
The single particle Hamiltonian can be diagonalizod nu- 
merically, and it is therefore possible to calculate O for 
quite large system sizes with an error limited only by ma- 
chine accuracy. To discuss the finite size behavior of the 
overlap and hence the OE, we calculated O and AE for 
non-interacting systems with M up to 1200. 

We expect O to be given by Eq. (30) with a j replaced 
by a. We can determine a by calculating O for differ- 
ent system sizes M[i], i = 1, 2, 3, . . ., taking the centered 
differences 



^Or. ^ _\n{0{M[i + 1])} - \n{0{M[z - 1])} 



\n{M[i + l]} -\n{M[i-l]} 



(35) 



and extrapolating by fitting a power law a + bx+cx^ with 
X = 1/M in some window of small 1/M. 

Alternatively we can calculate a using the results of 
BCFT Eq. (10). For non-interacting electrons BCFT 
predicts that the finite size corrections are given by in- 
teger powers of 1/M [24,25]. Thus we calculate AE for 
different M\i], take the centered differences 
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FIG. 1. The non-interacting orthogonality exponent ex- 
trapolated from a finite size scaling of the overlap and the en- 
ergy difference for a site and a hopping impurity with e/\t\ = 3 
and b = 0.1. 



a obtained from Eqs. (16) and (18), and the absolute 
extrapolation errors of a finite size scaling in different 
1/M windows are given in Table I. Fitting only the data 
with M[i] e [200, 1200] we obtain numerical values for 
a with a very high accuracy. For M[i] € [50, 200] which 
is comparable to the range accessible by the DMRG the 
relative extrapolation error is between 10~^ and 10~^, 
depending on the impurity type and whether the extrap- 
olation is of a'-' [i] or a'^^ [i] . Generally the accuracy is at 
least one order of magnitude better for the site impurity 
than it is for the hopping impurity. We therefore con- 
clude that while the hopping impurity has the advantage 
of generating only backward scattering, its finite size cor- 
rections are much more important than those for the site 
impurity. As there is no reason to believe that the ex- 
trapolation error for interacting electrons is smaller than 
that for non-interacting electrons this observation will be 
important in what follows. We also notice that the value 
of a obtained by the extrapolation of q;^^[z] is at least 
one order of magnitude better than the one from Q!'^[i]. 



1 

2ttv 



AE{M[i + 1]) - AE{M[i - 1]) 



(1/M[z + 1])-(1/M[z-1]) 



(36) 



and once again fit the data with a quadratic polynomial 
in 1/M. 

Fig. 1 shows a'~'[i] and a^^[i] for the site and hopping 
impurities in the strong impurity limit with £/\t\ — 3 and 
& = 0.1 at half filling. The results discussed in this sec- 
tion are quite general and not limited to the case of strong 
impurities. Prom Fig. 1 it is obvious that we can confirm 
the prediction of BCFT Eq. (10). The exact values of 



B. Accuracy of the DMRG 

The DMRG algorithm is a real space blocking method 
that works particularly well for the calculation of the 
ground state energies and wave functions of ID electronic 
lattice models. The accuracy of the method is principally 
limited by the number m of states that are kept after the 
diagonalization of two joint blocks. A compromise must 
be reached between keeping m small in order to reduce 
the numerical effort, and retaining enough states to give 
a good description of the system built out of the blocks. 
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FIG. 2. The overlap as a function of 1/m. The full symbols 
are the exact result for F = obtained by diagonalizing the 
one-particle Hamiltonian. For comparison, the full diamond 
is shifted by 0.014. The open symbols are the DMRG results. 
The circles are the overlap for M = 50, V — and a hopping 
impurity with 6 = 0.1. The squares are the data for M = 50, 
V/\t\ = 1 and a bond impurity with 6 = 0.1. The data are 
shifted by 0.03026. The diamonds and triangles arc the data 
for the same parameters but M = 70. The diamonds are 
shifted by 0.014 and the triangles by 0.0474. In the legend V 
is measured in units of \t\. 



obvious from Fig. 2 that the behavior of 0{M, m) is quite 
similar to that for V — 0. 

The same problem occurs for the other types of im- 
purity and in calculations of a from AE, as AE{M, m) 
always approaches A£'(M, oo) from below. The absolute 
error in AE is of the same order of magnitude as the one 
in O. 

The best way to overcome this difficulty would be an 

extrapolation from finite m, to m ~ oo. Unfortunately 
the functional dependence of 0{M,m) and AE{M,m) 
on m is not known. The numerical data for 0{M,m) 
and AE{M, m) show signs of a power law dependence 
on 1/m with an exponent of the order of 4.5, but we 
found it impossible to extract a convincing extrapolation 
law. By choosing very large m we instead ensure that the 
systematic error is always insignificantly small. To this 
end wc keep m = 600 states per block for M £ [6,48], 
m = 800 for M e [50,58], m = 1000 for M e [62,80], 
m = 1100 for M e [84,90], and m = 1200 for M = 100 
for all data sets shown in the next section. For these m 
the value of 0{M, m) - 0{M, oo) for M = 100 is smaller 
than 10~^, and it is even less for smaller values of M. As 
we will see in Sec. IV, this level of accuracy is needed to 
obtain meaningful results. 

As can be seen in Fig. 2 the error for the interacting 
data is larger than for the non-interacting ones. This 
observation holds for all M, AE and the other types of 
impurities. 



As the non interacting case is not a special situation 
for DMRG, a comparison with the "exact" results for 
O calculated as in Sec. Ill A gives us a reliable measure 
of its error. We therefore used the DMRG to numeri- 
cally calculate 0{M, m) for the site and hopping impu- 
rity with different impurity parameters, different M and 
m. For the same impurity parameter and M wc calcu- 
lated 0{M, oo) as in Sec. Ill A. As a generic example 
we show in Fig. 2 0{M, m) for a hopping impurity with 
V = and 6 = 0.1 as a function of 1/m for M = 50 and 
M = 70. Additionally the figure shows 0{M,m) for the 
bond impurity with V/\t\ = 1 and b = 0.1. These are the 
parameters used in Ref. [30]. To present the curves for 
the different M and V in one figure the data for V = 0, 
M = 70 and V/\t\ = 1. M = 50,70 are shifted as indi- 
cated in the caption. It is very important to notice that 
for all M the overlap approaches the m = oo limit from 
above and that 0{M, m)~0{M, oo) is increasing very fast 
with M. To calculate a from 0{M, m) we have to take 
a numerical derivative (see Eq. (35)). If we fix m and in- 
crease M, 0{M, m) is too large and 0{M, m) - 0{M, oo) 
increases exponentially with M. The numerical differen- 
tiation in Eq. (35) reinforces this effect and we recover 
values of a[i] which are systematically too small. For 
large M[i\ this leads to a significant error in a[i\. 

Although we do not know 0(M, oo) for V ^ it is 



IV. NUMERICAL RESULTS FOR INTERACTING 
ELECTRONS 

In this section we present and interpret numerical re- 
sults for the three chosen types of local impurities. In the 
discussion we will distinguish between the weak impurity 
limit where we expect to find results given by the Eqs. 

(29) and (32), the weak hopping limit where we calculate 
the overlap with the OBC ground state and expect Eq. 
(34) to hold, and the strong impurity limit. 

A. Weak impurities 

For weak impurities and attractive interaction pertur- 
bation theory predicts that the backward scattering con- 
tribution to the overlap tends to a constant (sec Eq. (32)). 
The remaining M dependence of O is given by the for- 
ward scattering contribution. To check this, and Eq. (29) 
for aj{V), we calculated O and AE for several values of 
y < and weak site impurities. From these data we 
fomid a'-''[i] and Q!^^[i] as determined by Eqs. (35) and 
(36). Fig. 3 shows the results for the parameter sets 
V/\t\ = -0.75, e/\t\ = 0.05, V/\t\ = -1, e/\t\ = 0.25 and 
V/\t\ = —1.5, e/\t\ — 0.1. Assuming that at large M the 
finite size scaling is given by integer powers of 1/M Eq. 

(30) and not an anomalous scaling of the form Eq. (32), 
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FIG. 3. The orthogonality exponent for weak site impu- 
rities and attractive interaction. The stars are a*^ and the 
circles axe a^^ for Vl\t\ = —0.75, e/\t\ — 0.05, the crosses 
are a° and the squares are a^^ for V/\t\ = —1, e/|t| = 0.25 
and the plus signs are a'^ and the diamonds arc for 
V/\t\ = -1.5, s/\t\ = 0.1. The full symbols represent the 
related af. In the legend V and e are measured in units of 
\t\. 



we can extrapolate the data by fitting polynomials up 
to second order in 1/M. As the fitting works very well 
this procedure is justified. In Table II we list a'^^'^'^^ (0) , 
a^{0), a^{V), and the extrapolated a{V) for a fit with 
1/M G [0.01,0.03]. Taking into account the extrapola- 
tion error and the fact that the Born approximation is 
even in the non-interacting case only an approximation 
(see Table II), we conclude that our data are consistent 
with the predictions of perturbation theory and the RG. 
There is no backward scattering contribution to a and 
af{V) is given by af{V). We do not believe that the 
difference between the analytic values of aj{V) and the 
extrapolated a can be explained by a residual backward 
scattering contribution. The difference is much smaller 
than the backward scattering contribution for V = [as 
discussed following Eq. (16) at{V = 0) = a{V = 0)/2] 
and for some of the parameter sets the extrapolated a is 
smaller than the analytical {V). As a'~'[i] and a^'^[i] 
tend to the same limit as M oo the data are consistent 
with the results of BCFT (Eq. (10)). 

We would like to emphasize that for weak impurities 
it is important to calculate O and AE with a very high 
numerical accuracy. The difference between O and 1 is 
very small, e. g. for V/\t\ = —0.75 and e/\t\ = 0.05 this 
difference is of the order of lO"**. Therefore results are 
only meaningful if O and AE have a much smaller error. 
As discussed in the last section the error in O and AE 
in our calculations is less than 10~^. 

The next goal is to confirm Eq. (32) for weak hop- 



FIG. 4. Scaling behavior of the overlap for weak impurities 
and repulsive interaction. The symbols are the DMRG results 
for parameters as indicated in the legend. The solid lines are 
power laws with exponent 2 — 2K = 0.25 and the dotted lines 
have the exponent 2 — 2K = 0.5. In the legend V and e are 
measured in units of \t\. 



ping impurities, i. e. small 1 — 6, with attractive and 
repulsive interaction and weak site impurities with re- 
pulsive interaction. In the case of the hopping impu- 
rity there is only backward scattering, and O is given by 
Ob- For site impurities we have to distinguish between 
the conventional weak forward scattering contribution to 
O ~ —af{V) In (M) and the anomalous scaling given by 
Eq. (32). We have already confirmed that the analytical 
result Eq. (29) for aJ{V) is a good approximation for 
small £/\t\. Therefore we can extract the anomalous be- 
havior by subtracting —aJ{V)ln{M). In the following 
O always stands for the anomalous contribution. 

As the constant ai{V) in Eq. (32) is unknown we ana- 
lyze O by calculating centered differences of O and expect 
from Eq. (32) 

^ 0(M [» + 1]) - 0(M[i - 1]) 
ln{M[i + l]}-ln{M[i-l]} 

= a2{V)-^{2-2K)M^-^^ + ... . (37) 

Vp 

In lowest order in perturbation theory Wb/vp for the site 
impurity is given by e/vp and for the hopping impurity 
by 1—6. As discussed in the Appendix 02(1^) (2— 2/f) > 0. 

For repulsive interaction {K < \) we expect x to in- 
crease as a power law. This power law holds for M with 
{Wb/vpYM'^-'^'^ < 1 but M > 1. For larger M the ex- 
pansion Eq. (32) breaks down and for smaller M the cor- 
rections symbolized by the dots become important. Fig. 
4 shows a double logarithmic plot of x as a function of M 
for weak hopping and site impurities with different 6, e 
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FIG. 5. Scaling behavior of the overlap for weak hop- 
ping impurities and attractive interaction. The symbols are 
the DMRG results for parameters as indicated in the leg- 
end. The solid line presents a power law with exponent 
2 — 2K = —0.3834 and the dotted line has the exponent 
2 — 2K = —0.7100. In the legend V is measured in units 
of \t\. 



and V. The lines are power laws with the expected expo- 
nent 2 - 2K. For V = 0.445 Eq. (13) gives 2-2K = 0.25 
and for V = 1, 2 - 2K = 0.5. For large M the DMRG 
results approach the predicted power law behavior. As 
we have chosen small 1 — 6 and e/\t\, the power law holds 
for all accessible system sizes with M > 40. 

For attractive interaction {K > 1) and weak hop- 
ping impurities Eq. (37) predicts that x tends to zero 
as M^~^^. This power law behavior should hold as long 
as2-2K > -1, i. e. for 1< < 3/2. For larger K the 
1 /M correction included in the terms symbolized by the 
dots dominates the large M behavior. 

If we calculate x/i^b/vp)^ the right hand side (rhs) 
of Eq. (37) is independent of (Wh/vp)"^, i. e. we expect 
the data for large M, fixed V but different Wb to collapse 
onto one curve. The lowest order expression of Wb/vp 
in terms of the microscopic impurity parameter is given 
by 1 - 6. We calculated x/(l - for V/\t\ = -0.5 
and V/\t\ = -0.8 each for b = 0.95 and b = 0.9. From 
Eq. (13) we obtain for V/\t\ = -0.5, 2 - 2^ = -0.3834 
and for V/\t\ = -0.8, 2-2K = -0.7100. Even though 
1 — & is very small the data points for the two different 
h do not collapse onto one curve and we find a deviation 
between the power law behavior we expect from Eq. (37) 
and the behavior of the data at large M (see Fig. 5). This 
shows that higher order corrections are of importance. 
From a numerical point of view it is impossible to work 
with smaller 1 — 6, as the maximal difference between 
the overlap O and 1 for 6 = 0.95 is only 1.5 • 10~^. For 
smaller 1 — 6 this difference would be even smaller and 



we would get very close to the numerical accuracy of O. 

We have to distinguish between two different sources of 
higher order corrections. The first of these arises in the 
mapping from microscopic to field theoretical impurity 
parameters. As discussed in the context of Eq. (20), the 
higher order corrections to a^{V = 0) for the hopping 
impurity are of the order (1 — 6)^. If we include this 
correction and make the identification 



M/2 

-f = (l-6)2[l + (l-6)] 



(38) 



the different data sets are only weakly 6 dependent. Even 
if we include these corrections the results for different h 
do not approach one another for large M. Instead the 
difference between the 6 = 0.95 and 6 = 0.9 data sets 
slowly increases as M gets larger. This can be explained 
by higher order corrections in Eq. (32) (see below). Fig. 
5 shows a double logarithmic plot of x/[(l — 6)'^(2 — 6)] for 
the above parameters. We also calculated xl^l with the 
exact non-interacting backward scattering phase shift 5b 
given by Eq. (18), but this gives no further improvement 
compared to the results shown in Fig. 5. We therefore 
conclude that even for relatively small 1 — 6 it is necessary 
to take higher order corrections in the mapping of the mi- 
croscopic impurity parameters onto the field theoretical 
parameters into account. 

The slight deviation of the numerical data from the 
expected asymptotic power law behavior (shown by the 
lines in Fig. 5) suggests that higher order corrections are 
also important in the field theoretical calculation lead- 
ing to Eq. (32). Unfortunately it is difficult to calculate 
these corrections analytically, even within the TL model. 
As, among the cases considered here, these effects seem 
to be important only for weak hopping impurities with 
attractive interaction, we will not pursue this question 
further in the present article. 

We also performed DMRG calculations for weak bond 
impurities. The results are equivalent to the results we 
obtained for hopping impurities. 

In summary, we conclude that our highly accurate 
DMRG data show the expected weak impurity behav- 
ior in most of the cases discussed. For site impurities 
and attractive interaction we have numerically verified 
the Eq. (29) for a^{V). The backward scattering con- 
tribution to a does vanish. For weak site, hopping, and 
bond impurities with repulsive interaction the data are in 
agreement with the scaling law Eq. (32). Only for weak 
hopping and bond impurities and attractive interaction 
we found small deviations between the scaling law Eq. 
(32) and the numerical results. In this case higher order 
corrections are of importance. 



B. Weak hopping in an open chain 

To verify Eq. (34) for our microscopic lattice model we 
calculated the overlap 0° between the chains with OBC 
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FIG. 6. Scaling behavior of the overlap 0° between the 
open chain and a chain with a weak hopping for repulsive in- 
teraction. The symbols are the DMRG results for parameters 
as indicated in the legend. The solid line is a power law with 
the expected exponent 2 — 2/K = —0.3217 and the dotted line 
has exponent 2 — 2/K = —2/3. In the legend V is measured 
in units of \t\. 



(hopping impurity with 6 = 0) and with a strong hopping 
impurity, i. e. small tb, for repulsive interaction. There 
is no forward scattering contribution to the overlap and 
0° is given by alone. Again it is very important to 
have numerical data for 0° with a large accuracy as the 
difference between 0° and 1 is of the order of lO""*. 

The constant aiiV) in Eq. (34) is unknown and as 
above we analyze 0° by calculating centered differences. 
For these we expect 

^ Q°(M[^ + l])-Q°(M[^-l]) 
ln{M[i + l]}-ln{M[i-l]} 

= a2(F)|(2-2/i^)M2-2/^ + ... , (39) 

to hold. As in the weak impurity case we have {V) (2 — 
2/K) > 0. Provided that 2 - 2/K > -1 and < 1 
(i. e. 2/3 < K < 1) x°W decays as a power law in 1/M, 
with exponent 2 — 2/K. For smaller K the 1/M cor- 
rection to Eq. (39) (among the terms indicated by the 
dots), dies away more slowly than the M^~^/^ term, 
and dominates the behavior at large M. Furthermore 
if we divide both sides of Eq. (39) by {h/t)^ the rhs is 
independent of tb. For the weak hopping case {h/tY is 
given by 6^. Therefore Eq. (39) predicts that the data 
for large M, different 6, but the same V collapse onto the 
same curve. Fig. 6 shows a plot of x° /^^ ^is a function 
of M for V/\t\ = 0.5 and V/\t\ = 1, each for h = 0.05 
and b = 0.1. The analytical exponents 2 — 2/K for these 
V are —0.3217 and —2/3. The lines are power laws with 
the expected exponents. 



FIG. 7. Behavior of a [i] for a repulsive interaction in the 
limit of strong bond impurity. The symbols are the DMRG 

results for parameters as indicated in the legend. The dotted 
line is a = 1/16. In the legend V is measured in units of \t\. 



As predicted by perturbation theory, the data for dif- 
ferent b but the same V asymptotically collapse onto one 
curve. For large M these curves show power law behav- 
ior with the expected exponent 2 — 2/K. The numerical 
results are consistent with Eq. (34), i. e. the OE a° tends 
to zero and the OBC ground state and the groimd state 
with a weak hopping are contrary to the non-interacting 
case not orthogonal to each other. 



C. Strong impurities 

The main question that we will address in this section 
is whether the backward scattering contribution to the 
OE for repulsive interaction tends to 1/16 as M oo, i. 
e. the open chain interpretation is valid. To avoid com- 
plications due to the forward scattering contribution to a 
we concentrate on the hopping and bond impurities. Re- 
sults for strong site impurities will be discussed elsewhere 
[39]. 

In Ref. [30] Qin et al. also discuss the question of 
whether ab tends to 1/16 for a bond impurity with 
b — 0.1, V/\t\ = 1 and up to 48 sites by plotting a'^[i] (sec 
Eq. (35)). As described in Sec. Ill A the extrapolation 
error in a'^[i] for a bond impurity is much larger than 
that for a site impurity, i. e. the finite size corrections 
are much more important. Because of the anomalous cor- 
rections expected for interacting electrons (see Eq. (32)) 
we cannot extrapolate by fitting a polynomial in 1/M. 
It is therefore questionable whether we can extract the 
M — > oo limit by plotting the bare data for a'-' [i] . 

Fig. 7 shows a'^[i] for several V and b. The squares 
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indicate data calculated using the same parameters as 
used in Ref. [30], but for much larger system sizes. Pro- 
vided that the open chain interpretation of the RG and 
perturbative arguments is correct, all the curves should 
converge to 1/16 in the limit M ^ oo. Contrary to Qin 
et al. we find that it is impossible to confirm or refute this 
prediction by plotting a.^[i\. Within the range of system 
sizes accessible, none of the curves saturates at 1/16 — 
finite M corrections dominate the scaling for all values of 
h. We note that the errors of the data presented in Fig. 7 
are in all cases smaller than the symbol size. It is known 
from BCFT that a'^[i] approaches 1/16 for 5 = [24], 
but even this simple case is not obviously confirmed by 
Fig. 7. The behavior of the data for the hopping impurity 
is similar. 

We have tested various different kinds of extrapolation 
procedure (c. g. second order polynomials and polynomi- 
als combined with terms which include anomalous scal- 
ing) as fits to the behavior of [i] in the limit M[i\ oo, 
but none of these proved satisfactory for either the bond 
or the hopping impurities. For a'-'[i\ the corrections to 
the M — > oo result are too large. In the following we 
therefore concentrate on Ai?. This means that we from 
now on assume that we can use the BCFT result Eq. (10) 
to calculate the OE a. In addition we will continue the 
analysis focusing on the more physical hopping impurity. 

At the RG fixed point 6 = 6* = the finite M correc- 
tions to A.E are given by integer powers of 1/M [24,31] 



A^(y,6*) = co + ci(V,6*)M-i 
+C2(V,6*)M-2 + , 



(40) 



with constants Cj, and the quadratic extrapolation of 

oo can be applied (see Eq. (36)). 



a^-^[z] for M[{ 
We extrapolated a 



for the hopping impurity as in 



the Figs. 1 and 3 and obtained afJ^^{V/\t\ = 0.5) = 
0.06252, atUV/\t\ = 1) = 0.06261, and atl{VI\t\ = 
1.5) = 0.06285. The absolute extrapolation error for 
non-interacting electrons and a strong hopping impurity 
(discussed in Sec. Ill A) is of the order of 10~^ (see Ta- 
ble I). The extrapolated a for V/\t\ = 0.5 is consistent 
with the value a = 1/16, within this error. For V/\t\ ~ 1 
and V/\t\ = 1.5 the difference between the extrapolated 
a and 1/16 is slightly larger. Taking into account the 
error due to the approximate nature of the DMRG and 
the fact that the extrapolation error might be larger for 
interacting electrons, we can still conclude that the re- 
sults for hopping impurities with 6 = are consistent 
with a = 1/16. The same holds for bond impurities. 

For 6 7^ we expect anomalous corrections in the finite 
size behavior of AE{V,b) and therefore to a^-^[i]. Un- 
fortunately we cannot extract the anomalous scaling be- 
havior directly from our data for a^^[i], as it is masked 
by the normal behavior Eq. (40). From the RG fiow 
equations given in Ref. [10], we expect that close to the 
b* = fixed point the effective b is given by bM~^, with 
—rj = 1 — 1/K, i. e. b fiows to zero as we increase the 
system size. The idea is now to replace b* in Eq. (40) by 



bM~'^ and expand the difference between Ai?(y, h) and 
Ai?(y, 6*) in bM~^. Before doing this it is important 
to notice that the constant co in Eq. (40) is a functional 
of the Fourier transform of the impurity potential Wb{k) 
and not only a function of the effective backward scatter- 
ing strength. Therefore we cannot introduce the effective 
b in this term and the leading term in the difference is 
an unknown constant cq. The expansion gives 



AE{V, b) - AE{V, b*)=co + b 



dc,{V,b) 



db 



6=6* 



+b 



dc2{V,b) 



db 



M-^-^-h... 



6=6* 



, 1.2 d^c,{V,b) 
^2 



(41) 



6=6* 



Taking the derivative (numerically the centered differ- 
ences) of this expression with respect to 1/M we arrive 
at 



AE{V,b) 



dci{V,b) 



+6 



db 
5c2(y,6) 



db 



d{AE{V,b) - AE{V,b*)} 
d{l/M) 

(l-h77)Mi-i/^ 
{2 + r])M-'^/^ + ... 



b=b' 



b=b* 



+ 2^ 



1 2 d-'c,{V,b) 



dip 



(1 + 27?)M2-2/^ + . . . 



(42) 



6=6* 



Provided we can use the scaling law for the effective 6 in 
the derivation of Eq. (41) we expect that our data can 
be described by Eq. (42). 

Fig. 8 shows —AE/{b\t\) as a function of M for a hop- 
ping impurity with V/\t\ = 0.5, V/\t\ = 1 and V/\t\ = 1.5, 
for 6 = 0.025, 6 = 0.05 and 6 = 0.1 (for a discussion of 
the data for 6 = 0.3 and 6 = 0.4 see below). The data for 
the different V and large M collapse onto a single curve. 
As we expect from Eq. (42), the 6 dependent corrections 
[e. g. the term in the fourth line of Eq. (42)] arc more 
important for weak interaction, i. e. K closer to 1, than 
for strong interaction. It is obvious from Fig. 8 that, for 
the accessible system sizes, the numerical data cannot be 
described by the leading term of Eq. (42) alone. Taking 
only the powers of M into account the term proportional 
to 6^ in Eq. (42) seems to be a more important finite size 
correction than that proportional to bM~^/^ . But as 
the data for small 6 and larger V have already collapsed 
onto a single curve for M > 50, we conclude that the 
term proportional to 6^ is already small and neglect this 
contribution in what follows. To verify that the data fall 
away as predicted we fitted the points for 6 = 0.025 and 
M > 30 with a function of the form 



/(M) = diM^-^/-^ + d2M-^/^, 



(43) 



with fitting parameters di, d^, and K. The fits are shown 
as lines in Fig. 8. The exact values for K are K = 0.861, 
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FIG. 8. Large M behavior of —A.E/{h\t\) for repulsive in- 
teraction and a strong hopping impurity. For each interaction 
strength V/\t\ = 0.5, V/\t\ = 1, and V/\t\ = 1.5 the figure 
shows data for b = 0.025, b = 0.05, and b = 0.1. The top set 

of curves is for V/\t\ = 0.5, the set in the middle for V/\t\ = 1, 
and the bottom set for V/\t\ = 1.5. The lines are fits to the 
data (for details see the text). The diamonds and triangles 
arc data for h = 0.3 and b = 0.4 both for V/\t\ = 1. In the 
figure Ab is measured in units of \t\. 



for V/\t\ = 0.5, K = 3/4, for Vl\t\ = 1, and K = 0.649, 
for V/\t\ = 1.5. These numbers must be compared with 
the exponents K = 0.873, K = 0.774, and K = 0.706 ob- 
tained by fitting the data. For Vl\t\ = 0.5 and V/\t\ = 1 
the agreement between the fitted and the expected K is 
satisfactory. For large interaction the agreement is not 
that good. In this case higher order corrections propor- 
tional to b seem to be more important. If we fit the data 
in a smaller M window, i. c. only for the few largest sys- 
tem sizes available, the difference between the fitted and 
exact K is decreasing as it should be. 

Wc conclude that for M — ^ oo the M dependent part 
of AE{V, b) approaches the M dependent part of the 
fixed point energy difference AE{V, b* = 0) in the way 
expected from RG. In particular we find no sign of a 
term of order in AEiV,b). If a^^{V,b) would be 
different from a^^{V,b*), AsiV^b) would contain such 
a constant given by 2TTv[a^^{V, b) - a^^{V,b*)]. By 
plotting AE{V,b) as in Fig. 8 we are able to detect 
such a constant if it is up to 100 times smaller than 
2nvF[a^^{V = 0,6) - a^^{V = 0,6*)]. As the numer- 
ical results for a^^^{V,b*) are consistent with a = 1/16 
it follows that a = 1/16 holds also for small b away from 
the fixed point value provided that the use of BCFT pre- 
dictions is valid. Wc wish to emphasize that we can give 
quantitative numerical evidence for the predicted RG flow 
of the ground state energy difference for a fixed 6, and not 
only a qualitative picture of the flow by choosing different 



6 as was obtained in Ref. [31]. 

We performed the same kind of analysis for the bond 
impurity. The results are similar to those presented for 
the more physical hopping impurity and the same con- 
clusions can be drawn. 

High accuracy DMRG then provides indirect numeri- 
cal evidence that, for repulsive interaction, the backward 
scattering contribution to the OE for the overlap between 
the PEC ground state and the state with a local impurity 
is independent of both interaction and impurity strength, 
and has the universal value 1/16. This result is consistent 
with those obtained by approximative analytical calcula- 
tions in the effective low energy field theory based on the 
open chain interpretation and discussed in the introduc- 
tion. 

Even though we have studied chains with up to 100 
sites, very small 6 and obtain DMRG data with a very 
high accuracy, we can only give indirect evidence for this 
universal behavior because finite size corrections are very 
large. To show directly that a'~'[i] (and a^^[i]) tend to 
1/16 would require much larger system sizes. 

Fig. 8 also contains data for 6 = 0.3, 6 = 0.4 and 
V/\t\ = 1. For the accessible system sizes the data points 
for these larger values of 6 are further away from the 
asymptotic curve onto which the data for smaller 6 col- 
lapse. However the points do seem to converge with this 
curve for large M. One can therefore speculate that the 
backward scattering contribution to the OE is also 1/16 
for larger values of 6. This is consistent with the open 
chain interpretation of the RG results. 



V. SUMMARY 

In this article we presented the results of an extremely 
accurate DMRG study of the orthogonality catastrophe 
in a ID lattice model of correlated electrons. Wc con- 
sidered three different types of impurity and a range of 
values of the interaction. By making a detailed compar- 
ison with exact resiilts for non interacting electrons we 
were able to establish the accuracy of our DMRG data, 
and arrive at a method of extrapolating our finite size 
results to the thermodynamic limit. 

For weak impurities we found that, in most of the cases 
considered, the numerical data can be well described by 
the results of a analytical calculation in lowest order per- 
turbation theory in the impurity strength. The pertur- 
bation theory is performed within the effective low en- 
ergy continuum field theory for ID correlated electrons. 
Higher order corrections are important in the case of 
weak hopping and bond impurities with attractive in- 
teraction. We have confirmed numerically that a weak 
backward scattering component to the OE scales to zero 
for attractive interaction and that the forward scattering 
contribution is modified by the interaction. The data of 
the overlap also show an interaction dependent anoma- 
lous scaling with system size which was predicted by per- 
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turbation theory in the (backward scattering) impurity 

strength. 

In the case of a strong impurity we were able to confirm 
indirectly the RG prediction that the backward scattering 
contribution to the OE tends to 1/16 in the thermody- 
namic hmit. This value is universal and independent of 
both impurity strength and the size of the interaction. 
Contrary to the conclusion drawn in previous studies by 
Qin ct al. we found that direct evidence for this behavior 
cannot be extracted from data for the overlap at the ac- 
cessible system sizes, because the finite size corrections 
arc very large. However we obtained quantitative numer- 
ical evidence for the scaling of the ground state energy 
difference predicted by RG. 

The open chain interpretation of perturbative RG re- 
sults leads to the prediction that ab{V > 0) scales to 1/16 
in the presence of an impurity with an arbitrarily weak 
backward scattering component (e. g. a hopping or bond 
impurity with b > 0.9). By close analogy ah{V < 0) 
should then tend to even for very strong bare back- 
ward scattering (e. g. a hopping or bound impurity with 
b < 0.1). We found it impossible to confirm these pre- 
dictions quantitatively on the basis of DMRG data for 
systems with up to 100 sites. The finite size corrections 
are in all cases too large, and the scaling too slow. Much 
larger systems are needed to make a reliable numerical 
determination of these issues possible. The results are 
however in qualitative agreement with this extension of 
the RG results. 
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APPENDIX 

In this Appendix we will briefly describe how to obtain 

the perturbative expressions Eqs. (32) and (34) for the 
overlap [6] . By Tig we denote a general Hamiltonian and 
by W the perturbation. The overlap between the ground 
state \Eo) of Hq and \E^) oiHo + W is given by 



0'^\{Eo\E^)\' = l 



E 

njtO 



{Eo\Ei)\ 



(44) 



where jiJ^) are the excited eigenstates of Ho+W. Insert- 
ing the expression for \E^) in lowest order perturbation 
theory in W gives 



E 



|(.Eo|W|^„) 

{En — Eof 



(45) 



By introducing a 5 function and taking into account that 
the smallest energy difference En — Eq is proportional to 
1/L we end up with 

« 1 - / '-^{Eq \W5 {y - [Ho - Eq\) W\ Eq) (46) 

We will now specify Ho and W to be given by the 
Eqs. (21) and (31). In this case the correlation function 
occurring in Eq. (46) is related to the 2kF component of 
the density-density correlation function in the TL model. 
For the TL Hamiltonian it is known how to calculate 
correlation functions like the one occurring in Eq. (46) 
[23]. Evaluating the expression gives 



-C2{V) 



-2K 



477^ Vp 



2K 



(47) 



with positive constants q. An explicit expression for 
C2(U) can be given, but in the following we only need 
to know that C2(U = 0) = 1 [12]. Eq. (47) is equiv- 
alent to Eq. (32). In the limit isT 1, i. e. F 0, 
(i2-2if _ - 2K) tends to In (L) and we recover 

Eq. (33). 

Taking the dual Hamiltonian describing an open chain 
and a representation of the weak hopping in terms of the 

dual field theory [10] instead of Eqs. (21) and (31) we can 
perform a similar kind of calculation and end up with Eq. 
(34). 

In the non-interacting case we can furthermore write 
down higher order corrections. They sum up to the ex- 
pected Ob ~ exp{— a6ln(L)} behavior and finite size 
corrections that are integer powers oi 1/L. 
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exact error from O error from AE error from O error from AE 

M € [200, 1200] M € [200, 1200] M € [50, 200] M € [50, 200] 



site 4.89322299 ■ 10"^ 1 • 10"'' 1 • 10"^ 2 • lO"'^ 5 ■ 10"'' 

hopping 4.76437497 ■ 10~^ 5 • 10"^ 3 • 10"'^ 2 ■ 10"'' 1 • lO"''' 

TABLE I. The first column gives the exact value of the OE a in the non-interacting tight binding Hamiltonian for a site 
impurity with e/\t\ = 3 and a hopping impurity with b = 0.1. The other columns are the errors obtained by extrapolating the 
data of Fig. 1. The second and third columns are the errors for an extrapolation with M € [200, 1200] and the fourth and fifth 

the ones for M e [50,200]. 



^exact(O) a^(0) af{V) ag,tr(V^) atftriV) 

V = -0.75, e = 0.05 3.1650 • 10"^ 3.1663 • 10"^ 3.7789 • 10"^ 3.8152 • 10"^ 3.9189 • 10"^ 

V = -l, £ = 0.25 7.8342-10-'' 7.9157 ■ 10"* 1.4072-10-^ 1.3862-10-^ 1.3991-10-^ 

y = -1.5, e = 0.1 1.2644- 10""* 1.2665-10"* 7.4598-10"* 7.4128-10"* 7.4294-10"* 

TABLE II. The OE for weak site impurities. The first column gives the exact non-interacting a, the second column the 
non-interacting Born approximation, the third column the forward scattering contribution of the interacting Born approxima- 
tion, and the fourth and fifth columns the extrapolated values from Fig. 3. The extrapolation is done by a quadratic fit with 
1/M € [0.01, 0.03]. In the Table V and e are measured in units of \t\. 
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